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ABSTRACT 


Preliminary results of an experiment to determine the 
accuracy and the effects of varying amounts of environmental 
data at the NAVTORPSTA, Keyport, Washington, Dabob Bay 
Range are presented. A 75-KHz source was tracked at six 
depths at each of six horizontal ranges. The data were 
analyzed by the NAVTORPSTA NUTRACK III isovelocity-layer 
computer program and by an isogradient-layer program 
developed here, STUTRACK I. The results show that a low 
signal-to-noise ratio is not a problem even at long distances 
and that an isogradient analysis yields superior results 
with less environmental data than does the isovelocity 


method. 
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ioe LN@RODUCTION 


A. DABOB BAY FACILITY 

A facility for the three-dimensional tracking of 
surface ships, submarines, and torpedoes is located at 
the Dabob Bay facility of the Naval Torpedo Station 
(NAVTORPSTA), Keyport, Washington. Chosen because of its 
favorable oceanographic features and its proximity to the 
_NAVTORPSTA, Keyport, Dabob Bay has been used for torpedo 
testing since 1949. It consists of a 250-KHz "high 
frequency" installation no longer operating, and a 75-KHz 
"low frequency" primary tracking system enveloping nearly 
the entire bay providing for underwater tracking in a 
volume approximately 30,000 feet by 4,500 feet by 600 feet 
deep. 

Dabob Bay, a branch of Hood Canal, is in a deep (600 feet) 
depression adjacent to the Olympic Mountains. The Quilcene 
River flows into the northernmost part and the Dosewallips 
mie saromrne Seuthermmosy.. — the sides of ¥the bay axe 
precipitous and predominately rocky , white the bottom is 
mud. Annual precipitation averages about 51 inches and 
while average snowfall is about 15 inches, any one snowfall 
seldom remains longer than a few days. Winds generally blow 
along the length of Dabob Bay; i.e., North and south. 
Southerly winds are uSually accompanied by clear skies and 


occur most often in summer. Normal wind velocities are 
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5 to 15 miles per hour. Tidal levels in Dabob Bay range 
from about -5 to +15 feet. 

Water temperatures change markedly from season to season 
and even from day to day. The changes are most noticeable 
near the surface but may be observed even at greater depths. 
Salinity varies considerably with depth as fresher surface 
water from terrestrial runoff overlays more saline water. 
The variation is also quite seasonal, reaching a maximum in 
ChHeespraige.. The main contributor to sound veloci*#y changes 
is water temperature, and the resulting sound velocity 


profile is representative of the temperature gradient. 


B. MOTIVATION 

NAVTORPSTA, Keyport, has developed a computer progran, 
NUTRACK III, which divides the bay into laterally homogeneous, 
isovelocity layers and traces a sound ray from a three- 
dimensional hydrophone array mounted just above the bottom, 
back Ge ivsmsource, 

An experiment was undertaken to provide comparison of 
Gie*lhree=dimensional posation of a%source with respect to 
a receiver as determined by various ray-tracing techniques. 
From this comparison it was hoped that an optimum ray-trace 
procedure would be indicated and that a better understanding 


of the significance of environmental variations would result. 


ol 





i. SeROCE DURE 


A. DETAILED RANGE DESCRIPTION 

The basic underwater tracking system on the 75-KHz 
fange s* comprised Of three components: hydrophone array, 
computer system, and vehicular instrumentation. The ship- 
board transducer emits a 75-KHz acoustic pulse in synchronism 
with a master clock at the computer site. This pulse is 


detected by each of the four hydrophones (RoR »R, and co, 


y 
in the array and mixed with a local oscillator having a 
frequency unique to that hydrophone. The outputs of these 
oscillators are preamplified and fed to a multiplexer, which 
sums the signals and transmits the composite signal via an 
Mmiderwavetscavle to the Computer site, where it is processed 
to determine the in-water transit time to each hydrophone 
(TT, »T and Teg 
1. Hydrophone Array Configuration 

The hydrophones in the NAVTORPSTA Range are arranged 
imma sidOtkh=paselaine System in groups Of four, each one 
1ocated on each of four adjacent corners in an imaginary 
30-foot cube, thereby defining the orthogonal coordinate 
system in which the measurements were made. These hydrophone 
arrays used buoyant spheres which exert an upward force in 
excess of two thousand pounds to keep them upright. The 
array is not free to rotate about the Z axis, but may tilt 


as much as two degrees from the horizontal in strong currents. 
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This tilt iS measured by Servo-pendulum transducers whose 
specified resolution is 0.00833 degrees, accuracy 0.025 
degrees, and linearity 0.05 percent over the full scale. 

The Range consists of six arrays (Fig. 1) aligned 
along a line 450 feet east and parallel to the Range axis, 
spaced at 6000-foot intervals, with the 00 array being the 
farthest north. The Range axis is 191° 18! 14.1" , 

2. Computer System 

The computer system, a Scientific Data Systems (SDS) 
model 920, iS installed on Zelatched Point and consists of 
three main subsystems: the signal processing subsystem, 
tne data collection subsystem, and the computer subsystem. 
The signal processing subsystem, the link between the hydro- 
phones and the data collection subsystem, processes the 
multiplexed Signal received from theiarrays, discriminates 
against unwanted Signals, and determines the in-water transit 
time to each hydrophone. The data collection subsystem 
establishes the master clock timing by which all timing is 
computed. It also prepares all the array signals for the 
computer subsystem, which calculates the tracked object's 
Dosition, prints the tracking data, and records the data on 
magnetic tape. For real-time examination, the data are 


PuiOotueG On X—-Y pilotters. 


peer Bh LMEN TALS P ROCE DURE 
Peace oO ACOUSTIC SOUrCeMWwLUN AaDDrOoOximately ten pounds 
of ballast was suspended from the stern of the sound boat 


mm cee U0, sb>5O sand -00—-foot depthstat six 


as: 





horizontal ranges (stations) from 900 feet to 4700 feet 
from the 02 hydrophone (Fig. 2). ‘The first station was 
located approximately 900 feet north of hydrophone array 02, 
and 450 feet east of the Range axis. Successive stations 
were located parallel to the axis of the Range; the second 
being directly over the array, the third 600 feet past it oo 
(diametrically opposite the first station), and then at about 
1500-foot intervals to about 4500 feet, beyond which signal- 
to-noise problems made further examination impossible. The 
source emitted one 1.3 millesecond 75-KHz pulse per second 
for sixty to ninety seconds at each depth at each station. 
The position of the boat on the Range was found in the 
horizontal from the Range's Autotape system, a portable and 
commercial microwave system used exclusively for ship 
tracking that ylelds accuracies of #3 feet. The position 
of the acoustic source with respect to the Autotape antenna 
was determined from the length of cable let out, the 
horizontal distance from the antenna to the stern of the boat, 
and the boat's heading. Tidal data were recorded and used 
for correcting depth coordinates to a mean low-low water 
(O-foot tide) reference plane. The array is known to be 
585 feet below MLLW, so that the depth of the source below 
the surface and the tide measurement allows calculation of 
Pic evcrpicaisaistancesbetweenmeuie source and array. All 
measurements were made with respect to the array's coordinate 


System, although the’array's orientation and location on the 
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Range were known, making it possible to transform this data 
to a true position on the Range. The X axis of the O02 array 
was 25° 24! 36,0" T, 

A NAVTORPSTA digital velocimeter that has an accuracy to 
+1.0 ft/sec and resolution of 40.7 ft/sec measured the sound- 
velocity profile at each station. NAVTORPSTA instrumentation 
at Dabob Bay measured and recorded the acoustic travel times 
from the source to each of the four hydrophones in the array 
to six decimal places for every source depth at each station. 
the three-dimensional position of the source with respect to 
the array was calculated from these travel times by the 
NAVITORPSTA NUTRACK III computer program. 

The data were gathered on 21 March, 1974 during a neap 
tidal period and while the sea surface was smooth and giassy, 
although there had been ripples earlier in the morning. The 
day was bright and sunny with scattered, thin, high overcast, 
the temperature was 55°F, and no noticeable wind was present. 
The shiny stainless steel transducer disappeared from view 
when about eight feet deep, indicative of thick biological 
concentrations (gelbstaffen) in the water. Rudimentary drift 
Medsitecnenussol DOdG mOvVeMemu relative to Surface water were 
taken by throwing paper into the water and watching its 
relative motion with respect to the boat, which indicated 


a value of approximately 0.3 ft/sec. 


Abe, 





IIl. THEORETICAL DEVELOPMENT 


The computer program analysis includes corrections for 
the fact that acoustic rays are not straight lines in a 
non-homogeneous medium. Two simple ray-tracing techniques 
exist: the water mass divided into layers of isovelocity 
water or into layers with a constant velocity gradient. 
The isovelocity method is employed in the NAVTORPSTA 
NUTRACK III program; while employment of the isogradient 


technique is one of the purposes of this present experiment. 


A. NUTRACK I1il 
1. Raw Array Coordinates 

The geometry involved in the calculations of tne 
raw array coordinates is shown in Fig. 3. As described 
eariien, the four hydrophones are on four adjacent vertices 
separated by a distance d, along the edge of the cube. The 
Sie Leroi CiCmarroaymecOOrdinalés etSmom bae Centerros ~The 
Speen Lune enesorvhocenal ceerdinates parallel to its¥edge. 
Let the transit times of The acoustic pulse from the 
tracked source to the c hydrophone and the Y hydrophone 


be T. and Ty respecvyi ve lye then 


(ame = (RAWY + se + (RAWX + dye + (RAWZ + g)° (1) 
and 
(oT, )° = @yie - g)° + (RAWX + dy2 + (RAWZ + Gye (2) 
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Poti? en Cmords perenecmpe lL Ween it and ie expanding, 


and combining terms yields 


2 F 
am (@ 
RAWY = Bal oat les Le) : (3) 


In similar fashion, analogous expressions may be 


obtained for the X and Z coordinates: 


2 
- CF a 
SEND Soe Ou i) aot) <a) 
and 
0 
IRENE, oe NE (5) 
where 
c = the average sound velcoity, about 4860 ft/sec; 
dg the distance between hydrophone’, 30 ft; 
Toot yetyet, Sera vye le >plmnes ama Seconds. 


The travel times are determined from a clock whose 
period is 0.4x107° sec. (2.5 mHz), which is sufficiently 
high to measure precisely the arrival time of a 75-kHz 


Signal (about 8 counts per quarter period). <A correction 


factor must be applied to correct for the array receiver 


integration (pulse width descrimination). Let 
I = the integration time, 800 us 
To = the measured travel time to the Ee hydrophone; 
qT. = actual travel time, 

so that 


ey 





| 
tt 
a= 
== 
KH 
Ad 


and 


rH 
+ 
bad 

It 


ee Ses 1 ane 
eae 21) 


Woou Syl aw 
cm y 


The relations for the raw coordinates are now 


2 
Cc t 
aq (Tyt+TY-KV(TE-T) (6) 


RAWY 


Where K = 21 = LA Pale sec and Po and Ty are the measured 


times. 
Similarly 
2 
= Ppempo = it 
RAWX 74 (7 Ty K)(T? ee Ce) 
and 
at 
—_- “ f a te 1 
RAWZ Te oo ee te K) (Tt ay, ; (8) 


ce. Tilt Correction to Determine Corrected 
Array Coordinates 


MHe Wwaw coordinavues Previously computed are 
remerenced CO 7a non—norizontal, reference plane. The varray 
is instrumented to sense the ee the array axes make 
iPePmneomiOrL Zonta called xALLLY or rotation about the’ x 


ewe LLL. CrerolLartonesavout the X axis, respectively. 
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These are given from 


WELT (TILTX - 505) x .00833 degrees 


and 


SATU GEE (TILTY ~ 505) x .00833 degrees 
where TILTX and TILTY are sensor outputs and 505 is an 
arbitrarily determined sensor reading for a perfectly 
horizontal orlentation. 

The correction of RAW coordinates for array tilt 


is merely a matter of geometry. Figure 4 shows that 
CORX = [RAWX][COS(XTILT) ] = [RAWZ][SIN(XTILT)]. (9) 


Since XTILT is small, it may be assumed that COS(XTILT) = 1. 


The relationship then becomes 


CORX = RAWX ~ [RAWZ][SIN(XTILT)] . (10) 
SamMilariy 
CORN paeR AW Ye RAW Zs oinG eT Le) (ele) 


and 


CORZ = RAWZ + [RAWX][SIN(XTILT)] + [RAWY][SIN(YTILT) ] 
ae 


These relationships are valid for both positive and negative 
tilt as the SIN term also changes sign to maintain the 


relationship. 


I, 





3. Ray Path Refraction 

In a homogeneous body or water, the minimum time 
path between an output transducer and receiving hydrophone 
is a straight line. When the speed of propagation is 
spatially varying due to temperature and salinity variations, 
refraction occurs and the minimum time path is no longer a 
straight line. ‘Corrections must then be made to determine 
the true position of the source. 

From the classical laws of Fermat and Snell, a 
minimum-time path is computed, working back in time from the 
array center in the direction of the "apparent" position and 
tracing the array through each velocity layer. The process 
is continued upward in space and backward in time until the 
measured travel time is consumed. To determine ray-path 
oie Cl OM Ge MS OlmGmsOCed ls vrarene atmwOnec—-TOOUsIlmcrements 
and averaged over ten feet. Since the array depth is 
585 feet, there are 58 ten-foot layers and one five-foot 
layer. The vertical angle of entry of the ray at the 
array, computed from the tilt~corrected array coordinates 


(Herm 4 Ss 


ly CORZ* le (13) 


2 2 2 


fh a=-es lh 
= CORK" + CORY + CORZ 


mhe@entry angile ne computed from the CORR array 
Socrdinatestis assumed. to be equal to the entry “angle for 
the minimum time path. This assumption is good if the 


sourcem ise tar enough so that the wave fronts are plane over 
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the dimensions of the array, and the speed of sound is 
constant over the dimensions of the array. 


The travel time, T, to the origin is computed from 


rT =e (14) 
where 
5 = [CORX* + CORY® + coRZ° 31/2 <1 5) 
and 
Ro = [(CORX + 15)* + (CORY +15)° + (corz + 15)2]?/? ; 
| (16) 


where ne is the slant range to the C array hydrophone, 
and S is the slant range to the array center. 

It may be seen that the assumption of an average 
velocity, c, in computing RAWX, RAWY, RAWZ is not sensitive 
as’ the ratio of S to Ro is what is really important and any 
small error will tend to cancel out. 


The horizontal distance is 


fa le 


H = [CORX® + CORY*] (17) 
The first layer processed is between the deepest 


deipeimein the velocity profile that is shallower than the 


artaye onge Gone array itselit; 


Gey 


where A, is the array depth. 


d 


Ai 





The incremental slant range is 


AZ, 


4S. = Sin ~ (19) 


The incremental travel time in this layer, used to increment 


Pic MOUUSeCe Crave me Gime ie 1s 


AS, 
AT =— , (20) 
ec 
i 
where Cy is the average velocity in the layer. This 


incremental travel time is used to decrement the pulse 


travel time T. The horizontal distance for this time is 


AH, = AS 


1 COS A, (213 


dl 


The values AH and AZ are accumulated as the time clock is 
decremented layer by layer. 
Using Snell's law, the change in angle due to 


refraction upon entry into a new velocity layer is 


COS A Cc 
O 


COSCAME+IAA) ae om 
COS A, + AA Cy 


or 


COS 42 i cos A = (22) 


ae 





where 


1 O 

and 

‘SIN A, = (1 - cos® Ne ey 
The thickness of the layer is 

AZ = Zo - 2, ; (24) 

This process is repeated and the original travel 
time T is decremented with the incremental layer travel 
times until T is consumed. The ray-bent X, Y, and Z 
(RBNX, RBNY, RBNZ) are 
= CORX 

RBNX = ( =3>) © An , (25) 

rBNY = ( SRE) san, (26) 
and 

RBNZ = © AZ (27) 


This raypath migration process assumes a laterally 
homogeneous medium within the neighborhood of each array. 


indemendenw velocity propillessare mot taken at each array. 





This-~w td 1. ieroduce small discontinuities in track as the 
tracked vehicle is passed from one array to another. 
4. Error Estimation 
Returning to Eq. 4, and assuming equal errors, 
AT, = AT. = At in time measurement; the maximum expected 


error in RBNX (AX) is found from 


2 
= = o- _ 2 
RAWX + AX xq [(T, + At) (T oe) A] 
BB rp? _ 2) 1 op ibe as IE fe | (28) 
od C Xx d Cc me 
to be 
°F 
= Se + 
AX = [T,, T, JAt (29) 
At large ranges, T_ »~x Dy» SO 
a) 246 
AX = ae TiAt ‘ 
isyobe e-T, =—n oo, USO 
AX e 
a ee (=) ie (30) 
Similarly, 
2G B © 
>. 2 (5) NG (31) 
and 
AZ .. £ 9 
= 2 (5) Ata or Ge 
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Assuming the maximum error in time is one-fourth of the 


period of the 75-kHz signal, AT. = 3.3x107° Seic, §SO 


.~ AY, AZ. -3 
te * a iL oe We 


B. ISOGRADIENT PROGRAM "STUTRACK I" 
An isogradient computer program, STUTRACK I, was 
developed for comparison with NUTRACK IIl. 
1. Linear Gradient Theory 
Figure 6 shows a highly exaggerated speed of sound 
profile which has been approximated by two layers each 
possessing a constant gradient. In each layer, the speed 


of sound may be given as 
e(Z) = c, + GAZ | (33) 


where G is the gradient; G is negative in the upper layer 


and positive in the lower layer. Applying Snell's Law for A, 





= A= CONSTANT , (34) 
one may see that the raypath curves downward in the upper 
layer and upward in the lower layer. Let the radius of 
curvature for the arc between the angles Ay and Ay in the 


upper layer be R. Then 





: 


AZ = R(COS A. -COS A) = R COS A. (—————= (35) 
al a? ed Cy 
However, from kq. 33, 
Cy) ~ 5 = -GAZ 
and with Eq. 34, 
A ee 
R = ¢ COS Eh | (36) 


Thus for each raypath the value of Ris a constant 
in a linear gradient layer, and the path is an are of a 
circle. The center of the circle lies at the depth for 
which A would become 90° if the sound profile for the upper 
layer were to descend to greater depths. 

Assuming that the initial angle of depression of the 
raypath is A., from the PeOMe Loy mOlmthne om ani Wi Lame amcor 
the change in depth and horizontal range become 


i, P= (COS A, - COS A) (37) 


Ren 
Ce TOLORS: A, 
and 


Ap Bae (SINA- SINA,) , (38) 


oe 
GeCos A, 


where 4G is the increase in depth and r 1s the increase in 
tranee ei womethe reference’ position, and A ts thie’ angle of 


depression of the raypath at this range and depth. 
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When AZ is known, the new angle A may be found by 


a= cos”? tcos aj(S4441)1 . (39) 
O 


The time-of-flight along the raypath may be 


simplified to 





A 
es dA 
oo ot GOS A? (40) 
O 
which upon integration yields 
“i COS A_.,1 + SIN Ag 
At = g inl(y-p-srra@ ‘cos a?! (41) 


Cee LU RACK Construction 

STUTRACK I (App. A) is designed to allow the deter- 
mination of a minimum amount of environmental information 
necessary for agreement with the NUTRACK III results. The 
program determines the tilt-corrected coordinates and the 
tay beau@mentry anglesinco the array if the®same» manner as 
WER ACiwt Lie puueuses Homes9 to calculate the exit angle 4A 
from the layer, Eq. 40 to calculate the horizontal range 
traversed, and the time-of-flight from Eq. 41. The thickness 
of the layer, AZ, is determined on a point-to-point basis 


from the depths from which the velocity profile was established. 


ei 





The values Ar and AZ are accumulated and the original travel 


time T is decremented as in NUTRACK III until 


TG >a Te tite ; 


where 
TG = = At 3 
and 


e = .0005 sec. 


This value of e¢ ensures an accuracy of 2% AZ and i Ar 
VWilwieiietserichess if 1G 38 withinge of T, the process 
stops and the X, Y, and Z coordinates for this isogradient 


technique (XG,YG,ZG) are calculated as in NUTRACK III, 


xq = Sone ZAr , (42) 
YG = ae Gaia (43) 
H 
and 
zq= SR bar . (4h) 
When 


TG >a +t 


Ene program returns to the previous layer and a bracket- 


and-halving process begins in AZ until the proper AZ is 


oye 





Loum emote sty =saneee seriterion imnjtime. XG, YG, and ZG 
are then found by Eqs. 42,83, “anid” 34. 

STUTRACK I is also designed to permit the utilization 
Of thei velocity—profa lie data in varying incremental step 
sizes to allow a determination of the effect of a varying 
amount of environmental data. For example, if velocimeter 
data is gathered in two-foot steps, an incremental step size 
(INC) of 1 would automatically set AZ = 2 ft, and an INC of 
10 would set AZ = 20 ft. Consideration of the source depth 
must be made prior to Ghe.“nvocatiwon sol oLMTRACK I in Sits 
present form, as the setting of too large an INC can create 
y AZ greater than the array depth of 585 feet, resulting in 
an error message. 

STUTRACK I is also designed to analyze the data witn 
the same isovelocity technique as NUTRACK III, but in 
varying incremental step sizes for comparison with the 
isogradient method. These output coordinates are XI, YI, 
and ZI. 

The inputs into STUTRACK I are: velocity profile, 
mICremene sme, billy correct loOn ys Lime Of day in hours, 
minutes and seconds, Ts Ty Ts Ta: Roe hoNieeeand hiBhize 
The output is six histograms of the differences between the 
NUURACK LE. rayobemt coordinates, amdeymeir respective 
values as determined by STUTRACK I using both the isogradient 
and isovelocity methods. These differences and the STUTRACK Il 


values are also printed. 
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cee oA ERNMENT AL RESULTS 


A. VELOCITY PROFILE 

Figures 7 and 8 illustrate that the velocity profiles 
for station 1 (taken at 1115 hrs) and station 6 (taken at 
1315 hrs) are equivalent below 85 feet. (See Fig. 2 for 
station positions.) The data for station 6, represented by 
the "+"'s, shows the "afternoon effect" of solar heating. 
Both profiles increase in velocity with depth at a greater 
rate than would be caused by pressure alone, implying a 
Similarly increasing temperature profile. The colder 
surface water was fresher, however, due to the heavy spring 


rainfall and runoff from the snow-capped mountain peaks. 


B. AUTOTAPE DATA 

The Range's microwave Autotape system was used to track 
the sound boat. Its precision was much Less than expected 
so that comparison with Range data was not attempted at this 
time and will have to wait until the Autottape data is 


Subjected to extensive analysis. 


Ceo Ayo LovOrF NUTRACK Tit RESULTS 
1. One-Way Transmission Loss 
The one-way transmission loss (TL) is determined by 


an assumption of spherical spreading and the relationship 
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TL(dB) = 20 LOG R + aR : (45) 
where 


attenuation constant in sea water 
for 75-kHz¥sound, 


Wy) 
tt 


lb 


0.03 dB/meter. 


Figure 9 shows that at short ranges (station 1) the TL is 
dominated by the spherical spreading term. Beyond a "critical 
range" of approximately 1000 feet for 75-kHz sound, the 
losses due to absorption of sound in sea water will overpower 
the spreading and the transmission loss will increase 
rapidly at greater distances. The TL from station 6 to the 
hydrophone array is approximately 105 dB, but is only 57 dB 
Prom st ation 1. 
2. RBNX, RBNY, RBNZ and Slant Ranges 

For Stations 1-25 and 6-50 

Figures 10-13 show the 80 data for station 1-25 
plotted versus time. Figures 14-16 represent the same data 
Awe eS ee O=-leeouG Wits linear draft correction. bisures 
i723) ere similar curves fior the 66 data for station 6-50, 
eexceD@ thateardqrift correction enor aied Only to the@silane 
range curve. These two stations oe chosen for preliminary 
investigation since station 1 was closest to the array while 
station 6 was farthest away. The RBNX, RBNY, and Slant Range 
scales for each station are nearly equivalent, but RBNZ data 
are plotted in a 5/1 ratio, since station 6 was five times 


CneuSianecanlrte Om SLaylon eines cases, station 1 data 


le 





show good point-to-point continuity while station 6 data, 
except for Slant Range, are very-erratic. Microstructure 
was examined as the cause of this variation using the 
procedure of Medwin [2], but was determined to be capable 
of producing a maximum fluctuation of only +1 foot, or only 
a small fraction of the actual variation. 
3. Transformation to Spherical Coordinates 

A transformation from the orthogonal coordinates 

RBNX, RBNY, and RBNZ to a spherical system was achieved 


using the following relations, 


R = (RBNX° + RBNY° + RBNZ°)272 (46) 
25) Beye 
p= smn) (RBNZ) 
and 
1 P RBNS 
9 = TAN (RET) 


where @ is the angle of elevation looking from the array to 
the source, and © is an arbitrarily defined azimuthal angle. 
A graph of @ versus 0 is @ Scatter diagram of the projection 
of the data onto the surface of a sphere of radius R. 

Figures 22-29 are station 1-25 plots of $ and 0 as 
calculated from: 


a. Ray bent coordinates (Figs. 22-24), 


peaeehay bent cCOordinates, amplified, and with linear 
iig@vincorrect lOmmindes .25—2 {in 


GeCUR coordinates, Mmplified (Figs. 28 and 29). 


Be 





Figures 22-24 show a consistent, nearly linear Grit t., With 
good point-to-point continuity. Figures 26 and 29 are 
identical, as should be the case since NUTRACK III assumed 
a laterally homogeneous medium, thereby keeping the ratio 
Y/X constant. Figures 25 and 28 show Similar but not 
identical variations, Which was also expected since the COR 
coordinates assumed a vertically homogeneous medium and the 
raybent coordinates consider the inhomogeneities. 

Figures 30-32 are station 6-50 plots of @ vs. time, 
O vs. time, and $6 vs. 0, respectively, as calculated directly 
from the raybent coordinates, and plotted to the same scales 
as Figs. 25-27. The variations are approximately the same 
in magnitude as for the same plots at station 1-25. However 
the scatter diagram for station 1-25 transforms into a five- 
foot square patch on the surface of a Sphere of radius 90 
feet. At station 6-50 this uncertainty grows with the ratio 
of the radii to a patch twenty-five feet square. 

4. Error Determination 

An error determination was made from the data from 
Figs. 12, 14, 15, and 17-19 by subtracting each datum from 
the average value drift-corrected coordinate and dividing 
by R. Variations from the averages were found to be less 


6 with a standard deviation (o} less than TO 


1 aia Seale 
This latter is comparable with the theoretical error deter- 
ML Tro peeonmoimse Cede | Ll—A—4 son the. basis Of an error in 


determining the arrival time which is independent of range. 
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Dee ec OMMARTSON OF NUTRACK TI TO STUTRACK I 
Comparison of STUTRACK I data in both the isSogradient 
and isovelocity modes with the NUTRACK III raybent coordinates 


was achieved by the following relationships 


RBNX — XG (47) 


DXG 
DYG = RBNY - YG 
DZG = RBNZ - ZG 

and 
Dae BN (48) 
DYI = RBNY - YI 


DZI = RBNZ - ZI 


Svevlonmme—JGmdata (Pigs. 36-40) are consistent with 
station 1-25 data (Figs. 33-35), but much more dramatic 
Since it was taken at a much greater range. These graphs 
are plots of DX, DY, and DZ for both the isogradient and 
isovelocity cases as a function of incremental step size. 
They show that the STUTRACK I resultsS are essentially 
equivalent to the NUTRACK III data when INC = 1. However 
as INC iS increased, the Mode oe ies technique gives 
consistently poorer results while the isogradient method's 
results remain adequate as INC increases. In all cases the 
STUTRACK I coordinates followed the NUTRACK III raybent 
coordinates point-for-point, but were biased to one Side. 
A few seconds of computer time was also saved as a result 


© the ase@gradient technique at a greater INC than for the 
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comparable calculations via the isovelocity method since 
fewer layers had to be exdmined to“acnieve the same, or 
Bet ters) recom tas 

FPigure 39 is a histogram of DXG, INC = 1, at station 
6-50. Its Gaussian nature is characteristic of all the 


histograms of DX, DY, and DZ at every INC tested, at both 


Stacwons. 








V. CONCLUSIONS 


A. NUTRACK III DATA 

The large variations in RBNX, RBNY, and RBNZ compared to 
the smoothness of the Slant Range may be explained by an 
examination of Eqs. 3, 4, 5, and 15. The individual 
Seordinaves ere dependence upon therdiiterence of To and Ty» 


a Ore T,- Since this difference is small, a small error 


y? 
in time will make a large fractional error in the computed 
Seerolniace,. melnemuncertainvy in signal discrimination, 
addressed in section III-A-4 and estimated as sane & Scce 
pocomnces important, Causing significant random errors from 
pulse-to-pulse. This error becomes very large for any 
position for which the wave arrives at the two hydrophones 
at nearly the same time. For calculation of the Slant Range 


the arrival time is used alone so the fractional error is 


much less. 


B. COMPARISON OF STATION 1-25 TO 6-50 

Comparison of station 1-25 data to that from 6-50 shows 
the latter to have greater spatial variability. The trans- 
formation to spherical coordinates shows this difference is 
Mov Proschumen the aazular display. Tnis IS what would bé 
SDC Clecael@etne UCIMemellOmsmwere the Important Sources OF 
Val Lan llityee this conclusion as upheld by an experimental 
error determination at both stations that agrees with the 


Pnieorepiacal error Induced by the one-fourth period 
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uncertainty in the signal arrival times. Out to the 
greatest distances studied (over twice that used in normal 
Range operations), the precision in measuring arrival times 
is unchanged even though the signal-to-noise ratio decreases 
Pe eo > Ode Lberotoroe tes conc ludede that for 
Condens ML ware temo nose Dievaldinogs at these. imesof «this 
experiment any attempt to improve the signal-to-noise level 
by increasing the source level will be nonproductive and 

. Wasteful. Improved Range accuracy can only be obtained by 
improving the accuracy of the time measurements (above its 
present value of about one-quarter period of the 75-kHz 


Signal) or increasing the baseline of the hydrophone array. 


C. EFFECT OF ENVIRONMENTAL DATA 

The Weraphs Of DX, DY; and DZ for both stations show 
COmeelot VC ly seuiaumnOl Bene =COndLUlOms. prevailing ag pie waime 
of this experiment the isogradient technique of STUTRACK I 
provides better results with less environmental data than 
does NUTRACK III. This is consistent in each oe the plots 
but is explicit in DZ at station 6, where the isovelocity 
layer method results in 4a ite er owen in D4 as the 
incremental step size (INC) is increased, but little change 
in DZ in the isogradient technique. Conclusive in this is 
that much less environmental data need be collected if a 


isogradient program were used. 


If 





VI. COMMENTS -AND S20SMMENDATIONS 


This thesis treated the NUTRACK III data as a comparator 
for analysis for both modes of STUTRACK I, but no conclusive 
evidence has been PaEne Tents determine which is more 
accurate. The Range's Autotape system was to have determined 
this, but was so imprecise that a decision was unreachable 
without detailed analysis of the Autotape results. 

In all Z calculations at both stations the source depth 
plus Z was less than the array depth by three-to-ten feet. 

No hypothesis has been reached for the cause of the 
disappearance of these few feet. 

Further study of the data is needed to determine if the 
results from the intermediate stations are consistent with 
Che@econclusicns reached here. Thengethese conclusions! should 
be tested using sound-velocity profiles representing other 


Dabob Bay oceanographic conditions. Special consideration 


should be given to the apparent inaccuracy in 2Z. 
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RAWZ SIN(XTILT) 


The SIN(XTILT) is positive when the X transducer is 
above the horizontal plane passing through the C 
transducer, and negative when the X transducer is 


below the horizontal plane passing through C. 


FIGURE 4a) SILT CORRECTION 
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